Benchmark on neutron capture extracted from {d,p) reactions 
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Direct neutron capture reactions play an important role in nuclear astrophysics and applied 
physics. Since for most unstable short-lived nuclei it is not possible to measure the (n, 7) cross 
sections, (d,p) reactions have been used as an alternative indirect tool. We analyze simultaneously 
**Ca(d,p)*®Ca at deuteron energies 2, 13, 19 and 56 MeV and the thermal (n, 7) reaction at 25 meV. 
We include results for the ground state and the first excited state of *^Ca. From the low-energy (d, p) 
reaction, the neutron asymptotic normalization coefficient (ANC) is determined. Using this ANC, 
we extract the spectroscopic factor (SF) from the higher energy (d,p) data and the (n, 7) data. The 
SF obtained through the 56 MeV (d,p) data are less accurate but consistent with those from the 
thermal capture. We show that to have a similar dependence on the single particle parameters as 
in the (n,7), the (d,p) reaction should be measured at 30 MeV. 
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I. MOTIVATION 

Reaction rates of capture reactions are a crucial in- 
put to astrophysical network calculations. In particu- 
lar, neutron capture reactions, which play a pivotal role 
in the astrophysical r-process nucleosynthesis, have to 
be known for nuclei between the valley of /3-stability 
and the neutron-drip line (see P, 0, Q and references 
therein). Typically the neutron capture rates for the r- 
process have been estimated using the statistical Hauser- 
Feshbach model, although this may be unreliable away 
from stability as the level density is low. (n, 7) cross 
sections for short-lived unstable nuclei cannot be mea- 
sured experimentally and have to be taken from theory. 
However, the production of unstable nuclei close to the 
r-process path has become possible in the recent years, 
and neutron transfer experiments like (c?, p) on these nu- 
clei are becoming more and more feasible as beam in- 
tensity continues to rise. Experimental programs using 
such transfer reactions to derive (n, 7) rates (i.e. [1,[E@|) 
have shown promising results. Theoretical work is needed 
to place this indirect method on firmer grounds. Cap- 
ture reactions can occur directly or through a resonance. 
This work focuses on the direct capture only and ana- 
lyzes for the first time {d,p) and (^,7) reactions simul- 
taneously using the combination of spectroscopic factors 
(SF) and asymptotic normalization coefficients (ANC) 
0> S [1] ■ In contrast to c harg ed particle capture which 
is mostly peripheral (e.g. OO]), neutron capture re- 
actions may contain an important contribution from the 
nuclear interior and consequently be very sensitive to the 
spectroscopic factor of the final state 0, d, • One nu- 
cleon SFs were first introduced into nuclear physics in 
the context of the shell model, where a sequence of or- 
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bitals is generated by a mean-field. SFs provide a mea- 
sure of the occupancy of these orbitals • Even with 
modern residual interactions, shell model predictions ap- 
pear to overestimate SFs for well bound systems when 
compared to experimental values [13, [3| • The cause for 
this disagreement is still not well understood and work 
along these lines is ongoing. SFs on unstable nuclei are 
usually extracted through reactions such as transfer or 
knockout [Tsj . Transfer reactions have been traditionally 
the prime method of spectroscopy in nuclear physics (see 
e.g. the recent compilation of {d,p) reactions [3l). In 
the standard analysis, the experimental cross section is 
compared with the predictions from the distorted wave 
Born approximation (DWBA) [l3| or the adiabatic dis- 
torted wave approximation (ADWA) [l^, and the SF is 
extracted from the normalization. At present, there is a 
large interest in {d, p) reactions, be it as an indirect mea- 
sure of (n,7) rates, as a testing ground for many-body 
structure models, or directly associated with stockpile 
issues [Tgj . In the conventional approach, the extracted 
SFs suffer from large uncertainty due to the ambiguity 
in the bound state and optical potential parameters. As 
example we note the recently measured {d,p) reactions 
on ^^Ge and ^*Se nuclei in inverse kinematics, which have 
been used to determine the neutron SFs of ^^Ge and ^^Se 
above the closed neutron shell at = 50. These SFs 
were used to calculate the direct (n, 7) capture cross sec- 
tions on ^^Ge and ^''Se [lO]. However, the measured low- 
energy {d,p) reactions are peripheral; this leads to large 
uncertainties in the extracted SFs due to the ambiguity 
of the bound state potential parameters, which can devi- 
ate significantly from the standard ones for neutron-rich 
isotopes. For all these reasons, the methodology used 
when analysing (d,_p) data has been recently revisited 
and questioned 0, [g] . As DWBA predictions are often 
the base for the extraction of the phenomenological SFs, 
sources of uncertainty in DWBA calculations need to be 
under control. Once one has ensured the validity of the 



one-step approximation f2]|, one needs to assess the un- 
certainty in the optical potentials [2^ and the bound 
state potential parameters [1]. Along these hues, a com- 
bined method of extracting SFs from transfer reactions, 
using the asymptotic normalization coefficients (ANCs) 
determined independently, was suggested and tested for 
different nuclei 0, H, Hi- Introducing the ANC of the 
bound state in the formulation, one controls the normal- 
ization of the peripheral part of the reaction amplitude, 
generally a large contribution to the transfer cross sec- 
tion. This allows to significantly reduce the uncertainty 
in the choice of the bound state potential parameters and 
test the choice of the optical potential or the assump- 
tions in DWBA. These same ideas can also be applied 
to direct capture reactions, breakup, and (e,e'p). The 
ultimate goal of this work is to prove the principle that 
{d,p) reactions can indeed be used to extract (n,j) rate, 
through a systematic methodology. Along these lines, an 
important step to validate the use of {d,p) reactions as 
an indirect tool to extract (n, 7) cross sections for exotic 
nuclei is to check for consistency between SFs extracted 
from {d,p) reactions and those from (n,7) on stable nu- 
clei, when there is a large variety of data to constrain 
the problem. We therefore benchmark the (d, p) method 
on ''*Ca(n, 7)''^Ca (ground state and first excited state). 
This choice is based on the very high quality data for 
thermal neutron capture, the well known neutron scat- 
tering length, and the number of {d,p) data sets on ^^Ca 
which all together present a stringent test for the applied 
reaction model. The paper is organized as follows. In 
section II, we provide some theoretical background and 
present the procedure used for extracting SFs from either 
{d,p) or neutron capture. In section III we include the 
details of the calculations and the data used, and present 
the results of our test case. Finally, in section IV we draw 
our conclusions. 



II. THEORETICAL CONSIDERATIONS AND 
METHODOLOGY 

In both A{d,p)B and A{n, ^)B, the main nuclear struc- 
ture input is the overlap function between the final state 
and the initial state its norm being the SF. This 

many-body overlap function is usually approximated by 
a single particle state ifAn{r) such that the many-body 
effects are hidden into the normalization factor: 

lUr)=S'/'ipAnir). (1) 

Here, S is the neutron SF. We emphasize that Eq.([T|) 
is an approximation, since the radial dependence of the 
overlap function and the nuclcon bound state wave func- 
tion can differ: the overlap function, as a many body 
object, includes, in addition to mean-field effects, the 
effects of residual interactions, which affect the nuclear 
surface region, the relevant region for direct reactions 
[23} . The single-particle function ipAn{r) is the solu- 
tion of the Schrodinger equation with a central poten- 



tial, typically a Woods-Saxon potential of standard ge- 
ometry. Asymptotically, the many-body overlap func- 
tion and the single-particle function do have the same 
radial behaviour. Then, defining C as the asymptotic 
normalization coefficient of the overlap function, and b 
the single-particle asymptotic normalization coefficient, 
= Sb"^. It is clear that, if one knew b then know- 
ing the ANC C would provide directly the SF S. How- 
ever, this is hardly ever the case. Ambiguities in the 
single-particle parameters introduce ambiguities in the 
SF which are not well controlled 0,B[ll. The DWBA 
amplitude for the transfer reaction A{d,p)B is given by 

Af^P =< 'il;'f\An\AV\ippn^'i'^^ >. The transition oper- 
ator is written in post- form AV = Vpn + VpA — UpB with 
Vij the interaction potential between i and j and UpB 
the optical potential in the final-state. The distorted 
waves in the initial and final states are -i/'l"'''' and ip^j: \ 
and ippn is the deuteron wavefunction. Similarly, the re- 
action amplitude for the (n, 7) process in first order is 

given by M"^ =< ^An\0\(p['^}_ff >, where O is the well 
known electromagnetic transition operator, and ipscatt is 
the neutron incoming scattering wave, usually calculated 
with the same potential that generates the final neutron 
bound state fAn ■ The amplitude for (n, 7) is given in 
the usual first order approximation and the operator is 
taken in the long- wavelength limit. The A(d,p)B DWBA 
amplitude depends on the optical potentials in the initial 
and final states and the bound state potential, while the 
A{n, "f)B reaction depends only on the Va„ potential for 
the bound and scattering states. In either case, a phe- 
nomenological SF is typically extracted through normal- 
izing the cross section to the corresponding data. As in 
0, Q , and for illustration purposes only, one can split the 
total amplitude into an asymptotic part Mg^t = bM^xt, 
corresponding to r > i?jv where the single-particle func- 
tion already behaves as a Hankel function, and the re- 
maining interior part Mint ■ The cross section is schemat- 
ically given by 

a" (X 5 |Ar [6]|2 = 1^1/2 Mf„J6] + CM:J\ (2) 

where a stands for (d, p) or (n, 7) . Here we express the 
dependence on the single particle parameters (rg , a) by 
the single-particle ANC b — b{ro,a). Eq. ^ brings out 
important intuitive physics: the normalization of the ex- 
ternal part of the reaction amplitude is governed by the 
ANC while the internal part is determined by the SF. If 
a reaction is completely peripheral, it is possible to ex- 
tract the ANC from the normalization to the data with- 
out any single-particle ambiguity (i.e., no dependence on 
b). In general, there are contributions from both the in- 
ternal and external region and it becomes important to 
fix the external part independently for an accurate de- 
termination of the SF. This is true for transfer and cap- 
ture reactions. It is generally said that transfer reactions 
are surface peaked, however the internal/external relative 
contributions can change significantly with energy. For 
sub-Coulomb transfer reactions, one is only sensitive to 
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the asymptotic part of the neutron wavefunction, and the 
ANC can be extracted virtuaUy without theoretical un- 
certainties. In 0,0], the important reahzation was that, 
introducing this independent ANC into the formulation, 
one could then extract a SF from a transfer reaction at 
higher energy, reducing significantly the uncertainty from 
the single-particle parameters. Direct radiative capture 
(n, 7) can be mostly peripheral if there is a centrifugal 
barrier in the initial state (e.g. ^'*C(n,7)^^C [ISl)- How- 
ever if there is no barrier at all (s— wave neutron capture), 
the interior contribution is important. At typical stel- 
lar energies of 10-300 keV, s— wave and p— wave capture 
cross sections are often comparable, see e.g. [11], so both 
the ANC and the SF are important ingredients to calcu- 
late the (n, 7) cross section. In this work we will anal- 
yse, en par, (d, p) reactions at several beam energies and 
the corresponding capture reaction on "^^Ca. By vary- 
ing the single-particle radius within a reasonable interval 
and fixing the diffuseness, we generate a set of single- 
particle functions for a range of single-particle ANCs. 
For each of these, we calculate the {d,p) and (n.,7) cross 
sections. Normalizing the peak of the (d, p) angular dis- 
tribution to the data at each energy, and normalizing the 
thermal neutron capture total cross section, we calculate 
phenomenological SFs. We choose {d,p) at sub-Coulomb 
energies, where the reaction is known to be peripheral, 
and normalize the theoretical cross section to the back- 
ward peak in the data. We include (d, p) at higher ener- 
gies and normalize the theory to the first forward peak 
in the data. In this way we obtain a function Sf^p{h) 
from each deuteron energy considered. We also extract 
S'"Jp(6) from the thermal capture. These Sexp{b) func- 
tions obviously vary differently with the single-particle 
ANC b. If there is consistency in the formulations, they 
should all intersect for a realistic &o- This procedure is 
better illustrated by looking at the ANCs 
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If the reaction is completely peripheral, Cexp{b) is con- 
stant. If there is contribution from the interior, then 
Cexp{b) has a slope. This slope will be more pronounced 
the larger the interior contribution. Our test case is *^Ca, 
which is part of the r-process chart fl*! and believed to 
be composed of the double magic ''^Ca core and a neu- 
tron single particle state, either P3/2 for the "'^Ca ground 
state, or pi/2 for the first excited state in ^^Ca. Direct 
capture is dominant for ''^Ca(n, 7)^^Ca because of the 
small reaction Q- value and the low level density (see Fig. 
4 of [1^). These properties of the stable neutron-rich 
''^Ca are close to the expected properties of nuclei close 
to the r-process path. For "^^Ca, direct neutron capture 
occurs via s— wave [1^ [13, [H, [1^ and has a significant 
contribution from the nuclear interior. Therefore it is an 
ideal test case since the result should be very sensitive to 
the ANC, but also to the SF. Another advantage of this 
choice is that the scattering length for ''^Ca — n has been 
measured with sufficient accuracy a = 0.356 ± 0.088 fm 
[30| which allows to determine the initial scattering po- 
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FIG. 1; (Color online) Angular distributions for 

*®Ca(ci,p)'*^Ca compared to data populating the ground state: 
Ed = 13 MeV (red dotted line), Ed = 19 MeV (green dashed 
line) and Ed ~ 56 MeV (blue dot-dashed line) . Also shown are 
the angular distributions to the first excited state: Ed = 13 
MeV (black thin dotted line) and Ed = 19 MeV (black thin 
dashed line). 



tential Vau with very minor uncertainties [26[. Neutron 
radiative capture on "'^Ca [H, [13, [H, [1^ and transfer 
^*Ca(d,p)''^Ca [U, [H, [sl] have been accurately mea- 
sured at several beam energies, for both the ground state 
and the first excited state. Thermal (n, 7) data at 25 
meV [1^, [l?! were measured with 6% accuracy, and sub- 
Coulomb {d,p) at 1.992 MeV [31j] was measured with 10% 
accuracy. In this work we also use {d,p) data at 13 MeV 
and 19.3 MeV from [sl and at 56 MeV from [sl. 



III. RESULTS 

Neutron capture El cross sections are determined for 
the set of bound state wavefunctions corresponding to the 
single-particle radius within the range tq — 1.05 — 1.65 
fm and fixed diffuseness a = 0.65 fm, with the Woods- 
Saxon depth Vws adjusted to reproduce the correct bind- 
ing energy of the final states in ^^Ca. The depth Vu,s 
for the initial distorted wave is fixed to reproduce the 
n — ^^Ca scattering length [s^]. For the transfer calcu- 
lations, following the procedure in [^, we use the Perey- 
Perey global optical potential [13] to define the optical 
potential in the exit channel and construct the deuteron 
potential in the entrance channel using ADWA [3], for 
all but the sub-Coulomb reaction (for the 1.99 MeV {d,p) 
calculations, we use the optical potentials presented in 
[3H although results are not sensitive to this choice - 
e.g. a variation of the potential depth by 10% changes 
the reaction cross section by only about 1%). ADWA 
includes deuteron breakup and thus goes well beyond 
one-step DWBA. The n-^®Ca and p-^*Ca potentials, 
needed to construct the finite ran ge a diabatic deuteron 
potential [s^ , are also taken from |34l ] . ADWA is devel- 
oped for reactions where the remnant VpA — UpB can 
be neglected. We have checked that this is the case 
for ^^Ca(d,p). For Vnp we use the Reid-soft-core [1^ 
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FIG. 2: (Color online) 5,^^(6) from ^^Ca.{d,p)*^Csi{g.s.) at 
Ed = 1.99 MeV (green dots), Ed = 13 MeV (red squares), 
Ed = 19 MeV (purple diamonds), Ed = 30 MeV (open 
circles) and Ed ~ 56 MeV (open triangles), and from 
**Ca(n, 7)'*®Ca((;.s.) at 25 meV (blue triangles). Also shown 
are the experimental uncertainties in the {d,p) reaction at 
1.99 MeV (dashed lines), at 56 MeV (long-dashed lines) and 
the (n, 7) reaction (solid lines). 
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FIG. 3: (Color online) Same as Fig[2]but now referring to the 
ANC CLp(&). ANC and SF are related by Eq. Q. 



although a simple Gaussian provides the same results. 
For the deuteron bound state we use [s^ and for ^^Ca, 
the same set of bound state wavefunctions used for the 
(n, 7) calculations. All calculations are performed using 
FRESCO [371. The shape of the angular distribution 
at 56 MeV [33| , in particular in the forward angle region 
where it is most important, was not reproduced (this was 
also the finding in [IBl)- For that reason, instead of the 
Percy and Percy parameterization, we used the Koning 
and Delaroche [38] with which a significant improvement 
was obtained for the ground state. Calculated ADWA 
angular distributions arc compared to the data in FiglTJ 
Results for the "^^Ca ground state are shown in Fig[2]and 
FigOfor Sexp{b) and Cl^p{h) respectively. As expected, 
the sub-Coulomb (d,p) reaction is totally peripheral. At 
Ed = 1.99 MeV, the reaction is below the Coulomb bar- 
rier {Vc B ~ 5 MeV) both in the initial and final channels 
(Q-value is 2.924 MeV). C^xpi^) remains constant over 



the broad interval of b: Cl = 32.1 ± 3.2 fm"! for the 

P3/2 

49Ca(5.s.) ^ "^^Ca + n, see FigH Consequently the SF 
behaves as 1/6^ as illustrated in Figl^l The extracted 
ANC is insensitive to the optical and bound state poten- 
tial parameters, the error bar coming from the systematic 
error of the data alone. At Ed — 13 and 19 MeV, the 
((i,p) reaction turns out to be also dominantly periph- 
eral. It is reassuring that the C^^pib) extracted at these 
energies [U are consistent with the sub-Coulomb result, 
corroborating the ADWA at these energies. However, 
we should note that, contrary to the 2 MeV case, these 
cross sections are sensitive to the choice of the optical 
potentials. The lack of b dependence makes it impossi- 
ble to determine the SF from the 13 and 19 MeV data. 
Increasing the beam energy changes the picture. Fig[2] 
shows clearly that the 56 MeV data has a larger contri- 
bution from the interior and thus is more sensitive to the 
single particle parameters. Expectedly, as the single par- 
ticle ANC increases, so does the external contribution, 
which explains the flattening of the (6) curve from the 
56 MeV (d,p) data. The error band is represented by 
the long-dashed lines in FigslHandO The cross sections 
at this energy were measured with 10% accuracy but we 
have added in quadrature another 10% due to the depen- 
dence on the optical potentials. From the joint analysis 
of the 2 MeV data and the 56 MeV data, we can extract 
a spectroscopic factor Sexp = 0.55 ± 0.25. The results 
for the ANC extracted from '^^Ca(n, 7)^^Ca(/7.s.) show a 
strong b dependence, in perfect agreement with [39[, con- 
firming an important contribution from the nuclear inte- 
rior. From the overlapping regions in Fig|31 we obtain 
60 = 7.80 ± 1.2 fm-1/2 and Sexp = 0.53 ± 0.11, in agree- 
ment with the value obtained from the 56 MeV transfer 
data. This bo = 7.80 fm~^/^ corresponds to tq « 1.45 fm, 
a radius larger than the standard value. The determined 
SF is lower than one could expect from the independent 
particle shell model, confirming the reduction of the SF 
previously seen [lj| , and now the ambiguity of the bound 
state potential parameters has been eliminated. One may 
ask whether there are deuteron energies 20 < Ed < 56 
(MeV) for which the (d, p) cross sections show the same 
sensitivity to the nuclear interior as the corresponding 
(n, 7). We find that Ed = 30 MeV provides an excel- 
lent match. ADWA is expected to perform well in this 
energy region, as the adiabatic condition is satisfied (see 
[121 for a successful application on ^^C). We thus con- 
sider the {d,p) reaction at Ed = 30 MeV and apply an 
arbitrary normalization to the cross section (open circles 
in FigEland FigO. Clearly illustrated in FigH C^xpib) 
obtained from the {d,p) reaction at 30 MeV has a simi- 
lar b dependence as the C^xpib) from the (71,7) process. 
We emphasize that a {d,p) experiment on ''^Ca at this 
energy would be very useful; plans for such an experi- 
ment have been initiated by the present study |4^. We 
also study the (d,_p) and (n,7) reactions to the first ex- 
cited state in '*^Ca. Results for Sexp{b) and C^xpib) show 
the same pattern as for the ground state. Because of 
the relation between ANC and SF we show only C^xpib) 
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FIG. 4: (Color online) Same as Fig|3] but now referring to 
first excited state in *^Ca. 



in FiglH We obtain from sub-Coulomb {d^'p) reaction, 
Cl^^^ = 9.30 ± 0.93 fm-i. Reactions at 13 and 19 MeV 
remain peripheral and the 56 MeV distribution is not well 
described by our model. Again, the (d, p) at 30 MeV and 
the (n, 7) show a similar sensitivity to the interior as for 
the ground state transition. Our sub-Coulomb and 
(n, 7) joint analysis provides 60 = 3.63^Q'5g fm~^/^ and 



Sexp = 0.71_o;22- Again, this 60 corresponds to ro « 1.45 



fm, showing consistency in the geometry of the neutron 
P3/2 andpi/2 orbitals. The large radii ro give further evi- 
dence for a neutron skin in neutron-rich calcium isotopes 

El. 



IV. CONCLUSIONS 

Our goal was to provide a proof of principle that (d, p) 
reactions can indeed be used to extract (n, 7) rates, 
through a systematic methodology. As sub-Coulomb 
or near-Coulomb-barrier (d, p) reactions are peripheral, 
they provide an ANC virtually free from optical poten- 



tial ambiguities. Whenever the (n, 7) reaction is periph- 
eral, this is the only necessary bound state information 
needed for calculating the neutron capture cross section 
at low energy. If the (n, 7) is not peripheral, one needs 
to know both, the ANC and the SF. We have shown 
for ''^Ca that SFs obtained from thermal capture and 56 
MeV (d,p) are consistent. Our results suggest that 
data at energies around iS^ = 30 MeV has a similar inte- 
rior contribution as the (n, 7) and could thus be a better 
tool for extracting the (n, 7) cross sections. In order to 
confirm this, experiments at 30 MeV providing cross sec- 
tions with better than 10% accuracy, would be very use- 
ful. Any prediction of neutron capture cross sections on 
unstable nuclei will be complicated by the fact that the 
neutron scattering length will not be measurable in most 
cases. This leads to an additional uncertainty because 
the scattering potential is not as well-defined as for the 
presented example ^^Ca — n. A rough estimate (see e.g. 
Fig. 5 of Ref. [2^) shows that this additional uncertainty 
remains below a factor of two for realistic variations of 
the neutron scattering potential. A detailed study of un- 
certainties of neutron capture cross sections for unstable 
nuclei will be given in a forthcoming paper. 

From this work, we also find that s— wave (n, 7) reac- 
tions are actually very well suited for extracting a SF. 
For the (n, 7) process, the transition operator is well 
known, the ambiguity in the optical potential at very 
low energies is negligible when imposing the correct neu- 
tron scattering length, and the ambiguity of the bound 
state potential parameters is greatly reduced by fixing 
the asymptotic part of the bound state, say through sub- 
Coulomb reactions. Hence, the joint analysis of low en- 
ergy (d,p) and thermal s— wave (n, 7) provides a powerful 
method to test microscopic structure models. This work is 
supported by the NNSA-DOE Grants DE-FG02-93ER40773 
and DE-FG52-03NA00143 with a Rutgers subcontract DE- 
FG52-06NA26207, and the National Science Foundation, un- 
der grant PHY-0555893. 
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